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ABSTRACT 

We use the Ito stochastic calculus to give a simple derivation of the 
Lindblad form for the generator of a completely positive density matrix evo- 
lution, by specialization from the corresponding global form for a completely 
positive map. As a by-product, we obtain a generalized generator for a com- 
pletely positive stochastic density matrix evolution. 
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Completely positive maps, and dynamical semigroups constructed from them, play 
an important role in the theory of quantum dissipative systems [1]. In the global case, the 
general form of a completely positive map X — > T(X) has the simple and intuitive structure 
[2] 

T{X) = Y J A n XAi , (1) 

with the A n operators indexed by n. In the infinitesimal case, the most general generator of 
a completely positive density matrix evolution has the so-called "Lindblad form" [3] 

Cp = -i[H, p] + ^2[v n pvl - -pv ] n v n - ^vlv n p] , (2) 

nGiV 

with H self-adjoint. Although one normally expects a simple and evident correspondence 

between the global and infinitesimal forms of a transformation, the original proofs of Eq. (2) 

in Refs. [3] follow a less direct route. Recently, Peres [4] has made the interesting remark that 

the connection between Eq. (1) and the v n pv\ term in Eq. (2) can be heuristically understood 

by identifying v n with a rapidly fluctuating part of A n , of magnitude (dt)% . Much the same 

physical idea is exploited in the book of Parthasarathy [5], to give a derivation of Eq. (2) 

using stochastic calculus methods. Our aim in this note is to sharpen Peres' observation by 

using the ltd stochastic calculus to obtain Eq. (2) directly as an infinitesimal specialization 

of the global transformation of Eq. (1). As a by-product, we will obtain a generalized form 

for the generator of an infinitesimal completely positive stochastic density matrix evolution. 

We start by specializing Eq. (1) to operators A n of the form 

A n =d n + u n dt + v n dWt , 

(3) 

A{ =d n + u{dt + vldW t n . 
Here the d n are positive real numbers; no extra generality is achieved by taking the d n 's as 

complex, since by redefining the u n ,v n a phase in d n can be transformed into an overall 

3 



c-number phase factor in A n , which does not contribute to Eq. (1). The dWp are Ito 
stochastic differentials [6] which obey the algebra 

dW t m dW t n = c mn dt , dW t m dt = , (4) 

with c a real symmetric and positive semidefinite covariance matrix. By appropriately nor- 
malizing the operators v n , the diagonal matrix elements of c can always be made equal to 
unity, so that we have 

c nn = 1 , all n . (5) 

The definition of Eqs. (4) and (5) includes as special cases that in which the Ito differentials 
are all the same, dWJ 2 = dW t , all n (for which c mn = 1, all m,n), and that in which the 
Ito differentials dW™ are all independent ( for which c mn = 5 mn ). 

Consider now the completely positive density matrix transformation 

P^P + dp = T(p) . (6) 

Substituting Eq. (3) into Eqs. (1) and (6) and using Eq. (4), we get 
P + dp = ^(d n + u n dt + v n dW t n )p(d n + u{dt + vldW?) 

neN 

= ^d 2 nP +J2 d n(v nP + pvl)dW t n + ( P U^ + U P +J2 VnPvDdt , 



(7a) 



with 



U = ^ d n u n . (7b) 



Equating the coefficients of p on the left and right hand sides of Eq. (7a) gives the condition 



neN 



while for the change in p we get 



dp = d ^( v nP + pvl)dW t n + {pU ] + Up+J2 v n pvi)dt 



(8) 



neN neN 

Let us now take the expectation value of Eq. (8) with respect to the stochastic process; 
since the operators u n , v n have no dependence on the Ito stochastic differentials, and since 
E[pdW"] = in the Ito calculus [6], we get simply 



dE[p] = E[dp] = E[p]U^ + UE[p] + v n E[p]vl dt 



(9) 



neN 



Now let us impose the condition that the density matrix must always have trace 
unity, which implies that Trdp = Tr E[dp] = 0. From Eq. (9), we get by cyclic permutation 
under the trace 



= TrE[p] 



(10a) 



which can hold for general E[p] only if the coefficient operator is zero, which implies that 



(106) 



neN 



Thus the condition that the infinitesimal transformation preserve the trace of p determines 
the self-adjoint part of U, while the anti-self-adjoint part can be an arbitrary operator iH, 
with H self-adjoint, so that U has the form 



U = -iH -^Yl v " Vn 



:n) 



Substituting Eq. (11) into Eq. (9) then gives for the completely positive, trace preserving 
infinitesimal deterministic evolution of E[p] the expression 



dE[p\ 
dt 



= -i[H,E\p]] + Y, 



neN 



v n E[p}vl - \vlv n E[p\ - ]-E[p\vlv n 



(12) 



which is the Lindblad form. 

Returning now to the stochastic evolution of Eq. (8), substituting the above results 
and the condition Trrfp = 0, we get the additional condition 

TrpJ2dn(vn + vi)dW t n = , (13a) 

which can hold for general p only if the operator coefficient vanishes, so that 

J2 d n(Vn + vl)dW t n = . (13b) 
n£N 

Multiplying by dW™ we get the equations 

d n(v n + vl)c nm = , all m e N . (13c) 

nSiV 

The condition of Eq. (13c) can be rewritten by noting that since c nm is real symmetric and 
nonnegative, it is diagonalized by an orthogonal matrix O mr to yield nonnegative eigenvalues 



c r , 



£ c nm O mr = O nr c r , (14a) 



with no sum over r on the right. Thus Eq. (13c) is equivalent to 



dn{v n + vl)O nr c r = 0, all r £ N . (146) 

nGiV 

For those values of r for which c r > 0, we can factor out c r to give the restriction 

d n (v n + vi)O nr = , all r with c r > , (14c) 

while for those r for which c r = 0, there is no restriction. Defining diagonalized Ito differen- 
tials dZ\ by 

dZ r t = dW t n O nr , (15a) 



with the inversion 

dW t n = O nr dZ r t , (156) 

we find that 

dZ q t dZ r t = dW t n dW t m O nq O mr = c nm O nq O mr dt = S qr c r dt , (15c) 

n,m£N n,m£N 

and that the restriction of Eq. (13b) takes the form 

£ d n (v n + vl)O nr dZ r t =0 . (15d) 

Since Eq. (15c) implies that for those r for which c r = the corresponding dZ[ is idem- 
potent, and hence vanishes, the reduced set of restrictions given in Eq. (14c) suffices to 
guarantee the vanishing of Eq. (15d) and hence the satisfaction of the original condition of 
Eq. (13b). We conclude that the completely positive, trace preserving stochastic evolution 
of p corresponding to our construction of Eq. (12) is generated by 

dp =^2 d ^( v nP + pv n ) dw t + 

neN 

subject to the restrictions on the positive real numbers d n and the operators v n of Eqs. (7c) 
and (14c). 

The simplest case of Eq. (16) is that in which the sum over the index set N contains 
only one term, so that Eq. (7c) implies d± — 1 and Eq. (5) completely determines the 
covariance matrix to be c 11 = 1. Dropping the superfluous index n, we get 

dp = (vp + pv*)dW t + 

with the restriction of Eq. (14c) simplifying to 

v + v^ = , (18a) 
7 



-i[ H > P\ + z^Kpvl - \v ] n v n p - \pv ] n v n ] 



dt , (16) 



-i[H, p] + vpv^ — -V s vp — -pv^v 



dt 



(17) 



which implies that 



v = —iK 



with K self adjoint. Substituting Eq. (18b) back into Eq. (17), we get 



dp = -i[K,p]dW t + 
= -i[K,p]dW t + 
while from Eqs. (3), (7b), and (11) we see that 



- l [H lP \ + KpK- l -K 2 p- l -pK 2 



dt 



i[H,p]--[K, [K,p]] 



dt , 



A = l- iKdW t + (-iH - ^K 2 ^j dt = exp[-iHdt - iKdW t ] 



(186) 



(18c) 



(19) 



Thus, in the case when the index set N contains a single term, the transformation of Eq. (16) 
reduces to an infinitesimal stochastic unitary transformation. 
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